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Abstract 

We analyze a generalization of the sine-Gordon equation in laboratory coordinates 
on the half-line. Using the Fokas transform method for the analysis of initial- 
boundary value problems for integrable PDEs, we show that the solution u(x, t) can 
be constructed from the initial and boundary values via the solution of a 2 x 2-matrix 
Riemann-Hilbert problem. 
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1 Introduction 



We consider the following integrable generalization of the sine-Gordon equation: 

u X x - u tt = (l + v(d x + d t ) 2 ) wxu(x,t), x,teR, (1.1) 



where u(x, t) is a real-valued function and v £ K. is a parameter — note that ( 1.1 ) reduces 
to the sine-Gordon equation in laboratory coordinates when v — 0. In terms of the 
'light-cone' coordinates (£, 77) defined by 



t = Z-v, 



equation ( |1.1| takes the form 



u$ v = (1 + ^9|)sin(u), 



(1.2) 



(1.3) 



Equation (1.3) was derived using bi-Hamiltonian methods in|Fokas| (1995). It is related 



to the sine-Gordon equation in the same way that the Camassa-Holm equation is related 
to the KdV equation. 

In this paper, we will assume that v < and for simplicity set v = — 1. For this 
value of equation ( 1.3 ) appeared in Sakovich and Sakovich| ( |2007| , where it was shown 
to be related, via certain transformations, to an integrable equation which describes 
pseudospherical surfaces introduced in IRabelo (1989). The inverse scattering transform 



(1ST) formalism on the line for equation (1.3) with v 
and FokaslipfjlOl). 



T was implemented in Lenells 



A method for the analysis of initial-boundary value (IBV) problems for nonlinear 
integrable PDEs was announced in Fokas ( 1997 ) and subsequently developed further by 



T 




Figure 1 The half-line domain ft with respect to the laboratory and light-cone coordi- 
nates. 



Fokas (2008). Here, we use this method to study equation (1.1 1 in 



several authors, see 
the half-line domain 

ft = {(x,t)\0 < x < oo, < t < T}, 

where T < oo is a given final time, see Figure [T] Given initial values at x = and 
boundary values at t = such that the corresponding IBV problem for ( 1.1 ) in the domain 
ft has a solution u(x,t), we show that u(x,t) can be constructed via the solution of a 
2 x 2-matrix Riemann-Hilbert (RH) problem. The main notable features as compared with 
other similar applications of the methodology of Fokas (1997) are: (1) The formulation 
of the RH problem suggested by the above methodology depends, in addition to the 
variables (x,t), on a function p{x 1 1) which is unknown from the point of view of the inverse 
problem. In order to formulate a RH problem whose jump matrix involves only known 
quantities, we have to reparametrize the x and t variables. A similar situation occurs in 
the analysis of the half-line problem for the Camassa-Holm equation, where however only 
the a;- variable has to be reparametrized, see Boutet de Monvel and Shepelsky (20081. (2) 
Only certain combinations of u and its derivatives can be recovered from the RH problem. 
Therefore, in addition to solving the RH problem, the reconstruction of u(x,t) involves 
finding the solution of a nonlinear ODE. Following the ideas of Lenells and Fokas (2010) 
we show that thi s OD E can be reduced to a Ricatti equation. (3) The adopted Lax 
pair for equation (1.1) has singularities at A = oo and A = 0, where A G C = C U {oo} 
denotes the spectral parameter. In order to define eigcnfunctions which are bounded on 
the whole Ricmann A-sphere, we will use two different representations of the Lax pair. 
These representations are suitable for the definition of eigenfunctions which are bounded 
near A = oo and A = 0, respectively. 

The analogous problem for the sine-Gordon equation on the half-line (i.e. for the 
equation obtained from (1.1) by letting v = 0) was investigated in Fokas ( |2004 2008); 



Fokas and Its (1992 1. We emphasize that although there exists (when v < 0) a Liouville 



type transformation relating equation (1.1) to the sine-Gordon equation (sec Sakovich 



and Sakovich (2007); Lenells and Fokas (2010)), the half-line problems for these two 



equations are not equivalent, since the Liouville transformation transformation distorts 
the shape of the domain ft. 

In section[2]we introduce a Lax pair for equation ( |1.1[ ) and define bounded and analytic 
eigenfunctions which are suitable for the formulation of a RH problem. The jump matrix 
of this RH problem can be expressed in terms of certain spectral functions, which are 
introduced in section [3] Finally, the main result is stated in section [4] 
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2 Bounded and analytic eigenfunctions 



The Riemann-Hilbert formalism for integrating a nonlinear evolution equation is based 
on the construction of eigenfunctions of the associated Lax pair. These eigenfunctions 
are joined together to a bounded and sectionally analytic function on the Riemann sphere 
of the spectral parameter A € C = C U {00} . A Lax pair suitable for the construction 
of eigenfunctions which are bounded near A = 00 was derived in |Lenells and Fo kas 



(2010). For the problem on the line, this Lax pair representation alone was sufficient for 



the formulation of a RH problem, since the eigenfunctions could be constructed using 
only the x-part of the Lax pair. For the problem on the half-line, the construction 
of eigenfunctions involves also the t-part of the Lax pair, which has a singularity at 
A = 0. We will therefore introduce another representation of the Lax pair suitable for the 
construction of eigenfunctions which are bounded near A = 0. Then, according to the 



methodology of Fokas ( 1997 ) , we will define solutions of these Lax pair representations by 



integration from three different corners of the spatial domain £1. The eigenfunctions which 
are bounded near A = and A = 00 will be denoted by {/ij}? and {^j}^ respectively. 
Together the Hj's and the Qj's can be used to formulate a 2 x 2-matrix RH problem. 



2.1 Lax pair representations 

Let 

°1 = 1 n . a 2 = 



— i 

1 



^3 



1 

-1 



and define m(x,t) by 



m(x, t) = 1 + (u x (x, t) + u t (x, t)y 



(2.1) 



Applying the change of variables |1.2[) to the Lax pair of equation (1.3 1 derived in Lenells 



and Fokas (2010), we find the following Lax pair for equation (1.1) 



4> x + i (Xp x - ^) a 3 (f) = W+<j) 
4> t + i (Xp t + §y <t 3 = W-<f>, 



W± — =F ^cr 3 ± i 



(2.2) 

C U {00} is a spectral 
.cos(u) — (itt + u x ) sin(u) . (u t + u x ) cos(u) + sin(u) 



where <fr(x, t, A) is a 2 x 2-matrix valued eigenfunction, A G 
parameter, the functions W±(x,t,X) are defined by 



-cr 3 ± i- 



mX 



-<T2 (2.3) 



i{\ + cosu)(m(u t + u x )(l — cosu) — msin(ti) + 2u xx + 2u xt ) 

4(1 ± cos u)m 1 

and p(x,t) is a real-valued function such that 

p x = -(1 - cosu)^/TO, p t — -(1 + cosu)y/m. 



(2.4) 



The equations in (2.4) are compatible since equation ( |1.1| admits the conservation law 

((1 — cosu)y/m) t = ((1 + cosu)y/m) x . (2.5) 
We choose p(x, t) so that p(0, 0) = 0, i.e. 



p(x,t) 



(0,0) 



[(1 — cos u) \fmdx' + (1 + cos u)y/rndt' 



(2.6) 
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Despite the form of the denominator of the last term in (2.3), the functions W± do not 
have singularities at points where 1 ± cosu = 0. Indeed, using equation (1.1), the last 
term on the right-hand of (2.3) can be rewritten as 

m(u t + u x ) + u tt + 2u tx + u xx 

1(1 =F cosu) G i , 

4m 

and this expression is manifestly nonsingular. 

The functions W± have the following properties: 



• W±(x,t, A) — > 0, x -> oo, 

• ^OW)- ^ 0(1) 0(1/ X)) 

• tr(W±(ar,t,A)) = 0, 

• Wl(x,t,X) = —W±(x,t, A), 



A — > oo, 



where A' denotes the complex-conjugate transpose of a matrix A. The last two of these 
properties ensure that the eigenfunction (f>(x, t, A) can be normalized so that 



det(4>(x, t, A)) = 1, 4>\x,t,\) = -1 (a;,*,A). 



(2.7) 



The Lax pair (2.2) is convenient for the definition of eigenfunctions which are bounded 
near A = oo. In order to define eigenfunctions which are bounded near A = 0, we 
transform the Lax pair as follows. Let / denote the 2x2 identity matrix. The gauge 
transformation 

4>(x,t,\)=g(x,t)i>(x,t,\), (2.8) 

where 
g(x,t) 



\/2(u x + u t ) 



transforms (2.2) into 



((y/rh - 1) cos(u/2) + (u x + u t ) sin(u/2)) I (2.9) 
+ i (—(u x + u t ) cos(m/2) + (y/rn — 1) sin(u/2)J a\ 



(2-10) 



where Vj = Vj(x, t, A), j — 1,2, are defined by 

V\ — i sin 2 (M/2)(cos(w) — (u t + u x ) sinu)A<73 

— i sin 2 (u/2)((u t + u x ) cos(it) + sinii)A(72 

1 

— - (sin(u) — u t (l — cos it) + u x (l + cos it)) o\, 
1 

V2 = -(— 3 + 2 cos(u) + cos(2u) — (u t + u x )(2 sin(w) -I- sin(2u)))Ao- 3 

— i cos 2 (u/2)((w t + u x ) cos(u) + sinM)Acx2 



(2.11) 



+ - (sin(u) — u t (l — cosw) + u x (l + cosu)) u\. 
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(x,t) 



(x,t) 



(x,t) 



Figure 2 The contours of integration for the solutions [12, o,nd [13 of (2.15) 



The form (2.9) of g is motivated by the fact that it diagonalizes the terms of highest 



(2-12) 



order as A — > of the Lax pair \2.2\ and that it satisfies 

det(g(x,t)) = 1, gi{x,t)=g-\x,t). 



The relations (|2. 12 ) ensure that the gauge transformation (2.8) preserves the properties 
det(ip(x,t, A)) = 1, ipi(x,t,\) = ^ _1 (:r,t,A). 



in (2.71, i.e. 



(2.13) 

The function g(x, t) is nonsingular as u x + u t — > despite the form of the right-hand side 
of (pj|. In fact, 



g(x,t) 



sin(tt/2) — icos(u/2) 
-icos(u/2) sin(u/2) 



as u x + ut — > 0. 



The functions V\ and V% have the following properties: 
• Vj(x, t, A) — > 0, x — >• 00, j = 1,2, 



„, . u /0(A) 0(1)\ . . 

• tr(7jf(a?,*,A)) =0, i - 1,2, 

• VJ t (a;,t,A) = -^(x,t,A), j = l,2. 

2.2 Eigenfunctions bounded near A = 



1,2, 



In this subsection, we define solutions of (2.10) which are well-behaved near A = 0. 
Introducing an eigenfunction /i by 



wc find that the Lax pair (2.10) becomes 

( + i\[ a 3,fj] = VlfX, 

\ [Mt-i{\+ gx) [0-3, M] = 
This can be written in differential form as 

d ( e *(ft-< A +&>*)%) = W, 



(2.14) 



(2.15) 



(2.16) 



where 03 acts on a 2 x 2 matrix A by 03 A = [0-3, A], and the closed one-form W(x, t, A) 
is defined by 

W = e''(*-( A+ sx)*)*»(Vida: + V^di)^. (2.17) 
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D, 



Figure 3 The sets {Dj}f in the complex X-plane. 



We define three eigenf unctions of (2.161 by 



{ij(x, t,X) = I + 



(x,t) 



(2.18) 



where (xi,ti) — (0, T), (x 2 ,t 2 ) — (0,0), and (2:3,43) = (00, t). Since the one-form 
W is exact, the integral on the right-hand side of (2.18) is independent of the path of 
integration. We choose the particular contours shown in Figure [2] This choice implies 
the following relations on the contours: 



(xi,ti) 
(x 2 ,t 2 ) 
(x 3 ,t 3 ) 



(x,t) :x'-x<0, 
(x,t) :x'-x<0, 
(x, t) : x — x > 0, 



t' - t > 0, 
t' - t < 0, 
t' - t = 0. 



(2.19) 



Letting 



r> 1 = |Ae 

D 2 = i\€ 
D 3 = {\E 

d a = <: c 



Im — > and Im ( A 
8A 



Im — > and Im ( A 

8A 



Im — - < and Im ( A 
8A 



Im — — < and Im ( A 

8A 



1 

8A 
1 

8A 
1 

8A 
1 

8A 



< 

> 

< 

> 



(2.20) 



the second column vectors of jj,x, fi 2 , 113 are analytic for A € C such that A belongs to D3, 
D4, and Dili D 2 , respectively, see Figure [3} Moreover, away from A = 00 where the Lax 
pair is singular, they have continuous and bounded extensions to the closures of these 
sets. We will denote these vectors with the superscripts (3), (4), and (12) to indicate 
these boundedness properties. Similar conditions are valid for the first column vectors. 
We obtain 

Mi — y^i i J , A*2 — ( M2 > M2 ) > A*3 — ^3 , J ■ 
The fij's satisfy 

det(fij(x,t,\)) = 1, tXj(x,t,X)^ = ^(a;,i,A) _ \ j = 1,2,3. (2.21) 



G 



Moreover, as A — > 0, 

L^(x,t,X),^ 2 \x,t,\)] =7 + 0(A), A^O, AeD 2 , (2.22) 
(/4 M) (aj,t,A),A4 3) (a;,t,A)) = J + 0(A), A -> 0, A e D 3 . 

The /ij's are suitable for the formulation of a RH problem except that they have singu- 
larities at A = oo. Our strategy is therefore to cut out a neighborhood of A = oo and use 
the Lax pair (2.2 1 to define eigenfunctions which are bounded in this neighborhood. 

2.3 Eigenfunctions bounded near A = oo 



The form of the Lax pair (2.2) is convenient for the definition of eigenfunctions which 

(2.23) 



are well-behaved near A = oo. Introducing an eigenfunction $ by 

= $ e -»(*P+ i £r) ff 3 



we find that the Lax pair (2.2 1 becomes 

f $ x +i(Xp x -±)[(T 3 ,$] = W + $, 

This can be written in differential form as 

d(V (Ap+ ^ } * 3 $) = W^, 

where the closed one-form Woo(x, t, A) is defined by 



Wo 



(W+dx + W-dt)<&. 



We define two eigenfunctions $2 and $3 of (2.25) by 

f{x,t) 

$j(x,t,\)=I+ / e'^P^^+^^Wooix'^^X), 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



where (22,^2) = (0,0) and (xs,ts) = (00, t). The functions $2 and $3 are the analogs 
of /j,2 and /Z3 defined in (2.18); the analog of fix is not needed since we only consider 
a neighborhood of A = 00. Choosing the integration contours in Figure [2j the integral 
equations (2.27) defining $2 and $3 become 

$2(2, t,X)=I+ f e l ( A(p(0 '°^ (a; ' t))+ ^)* 3 (^$2)(0, f , X)dt< 
Jo 

° ( ! ' — \ 

e i(A(p<* .t)-f(^))-^J*3 (w+$3)( ^ )t)A)dx / > 

The second column of the integral equation for $2 involves the exponentials 

e 2i(A(p(0,t')-p(x,t))+^±£) and e 2i(\(p(x',t)- v (x,t))-^) 

where the functions p(-,t) and p(x, •) are nondecreasing. Define R > by 



$3(3, t, A) = /- 



(2.28) 



2 * - f 

R = SUp 7 ^rrrrrrrrrrrrrrrrr; . 

0<t',i<T4 ]*,(! + cos(u(0,r))i/m(0, r)dr 



(2.29) 
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We will henceforth assume that 

cos(u(0,t)) ^ -1, 



< t < T, 



(2.30) 



so that R is finite. The relations (2.19) together with the definition of R implies the 
following inequalities on the contour (£2,^2) —> (x,t): 

t' -t + x 



lm(x(p(0,t')-p(x,t)) + 
lm[X(p(x',t)-p(x,t)) 



8A 

x 1 — x 
8A 



>0 if \\\>R, ImA<0, 



> if Im A < 0. 



(2.31) 



The first of these inequalities is a consequence of the estimates 

t' -t + x 



lm(X(p{Q,t')-p(x,t)) + 



8A 



>Im(X)(p(0,t')-p(0,t) 
R 2 



t! -t 



8A 



>Im(A)(p(0,O-p(0,t)) 1 



IAI' 



Im A < 0, 



which hold on the contour (^2^2) ~~ * 
We define the sets D 5 and £> 6 by 



As 



I A G C|lmA > and |A| > i?| , 
I A G C|lmA < and |A| > R,X . 



(2.32) 



The inequalities (2.31) imply that [$2)2 is bounded and analytic for A G D§. Similar 
considerations apply to the other column vectors of $2 and $3, and we deduce that $2 
and $3 have the boundedness properties 



). *.= (*M 6) ) 



$2 = ( 

Moreover, the functions $2 and $3 satisfy 

det(®j(x,t,\)) = 1, $j(x,t,\y = &j(x, t, A) -1 



3=2,3. 



(2.33) 



3 Spectral functions 

We define two 2 x 2-matrix valued spectral functions s(A) and 5(A) by 
fn(x,t,X) ^ ^ 2 (x,t,X)e-<^ x+ ^'> t > 3 S(X). 



(3.1a) 
(3.1b) 



Evaluation of (3.1 ) at (x, t) = (0, 0) and (x, t) = (0, T) gives the following expressions for 
s(A) and S{X): 



s(A) = /z 3 (0,0,A), 5(A) = /ii(0,0,A) 



M 2 (0,T,A) 



(3.2) 



The function ^(x, 0, A) is defined by the integral equation obtained by setting t = in 
(2.18). This integral equation involves the function V±(x,0, A) defined in (2.11). Thus 
s(A) is defined in terms of the initial data u(x, 0) and Ut(x, 0) alone. Similarly, /ii(0, t, A) 



8 



is defined by the integral equation obtained by setting x = in ( 2.18[ ) which involves 
V2(0, t, A). Thus 5(A) is defined in terms of the boundary data u(0, t) and u x (0, t) alone. 
We infer from ( |2.21| that 

dets(A) = l, detS(A) = l, (3.3) 
and that there exist functions a(A), b(X), A(X), and B(X) such that 



^-6(A) a(A)y ^-B(A) A(A) y 

Dehning the sets D[ and D' 4 by (see Figure [5]) 

D[ = Di \ D 6 , D' A = D 4 \D 5 , (3.4) 

i.e. D[ and Z?4 denote the sets D\ and Z?4 with a neighborhood of A = oo removed, we 
can state the following result. 

Proposition 3.1 The spectral functions a(\) and b(X) have the following properties: 

(i) a(A) and b(X) are analytic for X € D\ U D 2 and continuous and bounded for X € 
D[UD 2 . 

(ii) a(A) = 1 + 0(A), 6(A) = 0(A), A -> 0, A e D 2 . 
(in) |a(A)| 2 + |6(A)| 2 = l, AeR. 

The spectral functions A(X) and B(X) have the following properties: 

(i) A(X) and B(X) are analytic for X € T)\ U D 3 and continuous and bounded for 
X G D[ U D 3 . IfT < oo, then A(X) and B(X) are defined and analytic in C\{0, oo}. 

(ii) A(X) = 1 + 0(A), B(A) = 0(A), X ->■ 0, X e DiU D 3 . 

(in) A{X)A{X) + B(X)B(X) — 1 for AeRU {|A| 2 = 1/8} f/orAeC ifT<oo). 



Proof. The properties denoted by (i) and (ii) follow from the discussion in subsection 2.2 



and the observation that the definition of /ii(0, t, A) implies that this function has the 
following enlarged domain of boundedness: 

m(o,t,x) = (^ 24) (o,i,A),/4 13) (o,;,A)) . (3.5) 

The properties denoted by (iii) follow from (3.3). □ 
We will also need the spectral function Soo(A) defined by 

$ 3 (x, t, A) = $ 2 (x, t, X)e-^P +t -^^ SoD (X). (3.6) 



Evaluation of (3.6| at (x,t) — (0,0) yields 

*oo(A) = * 3 (0,0,A). (3.7) 

Just like s(A), the function Soo(A) is defined in terms of the initial data u(x, 0) and 
u t (x,0). Moreover, satisfies 

det Soo (A) = l, (3.8) 

and can be written as 



,(A) 



(A) MA) 



(A) aoo(A) 

where aoo(A) and &oo(A) are complex- valued functions. 
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Proposition 3.2 The spectral functions a 00 {\) and b oa (X) have the following properties: 

(i) 000(A) and b co (X) are analytic in D 5 with continuous and bounded extensions to 
A e D 5 . 



(ii) | aoo (A)| 2 + |6oo(A)| 2 = 1, 



ImA = 0, XeD 5 



Proof. Property (i) follows from the discussion in subsection 2.3 Property (ii) follows 
from (Ol). □ 



4 The Riemann-Hilbert problem 



In this section we use the eigenfunctions and {"Sjjf to formulate a RH problem 

for a 2 x 2-matrix valued function with jump contour shown in Figure [5j We will first 
formulate a RH problem for a 2 x 2-matrix valued function M(x, t, A), whose form is 
suggested by the methodology of Fokas (1997). However, it turns out that the jump 



matrix for this RH problem depends on the function p(x, t) which occurs in the Lax pair 
(2.24). The function p{x,t) is unknown from the point of view of the inverse problem, 



and thus the solution is not yet effective. We will overcome this problem by introducing 
new variables (y, rf) and formulating a modified RH problem for a 2 x 2-matrix valued 
function M{y,i}, A), whose jump condition is given explicitly in terms of y, n, and A. 



The solution u(x,t) of (1.1) can be recovered in parametric from the asymptotics of 



M(y,rj, A). A similar reparametrization of the RH problem occurs also in the analysis 
of other equations such as the Camassa-Holm equation and equation (1.3), although in 



those cases only one of the variables (x, t) has to be reparametrized, cf. Boutet de Monvel 
and Shepelsky|(|2008|;|Lenells and Fokasl (120101). 



4.1 RH problem for M(x,t, A) 

We seek a bounded and sectionally analytic 2 x 2-matrix valued function M(x, t, A), which 
satisfies a jump condition of the form 



M- (x, t, A) = M + (x,t,X)J(x i t,X), 



A e D + n D- , 



M 



M+, 
M-, 



A e D+, 
A e 



(4.1) 
(4.2) 



where J(x, t, A) is a 2 x 2-matrix valued 'jump matrix' and 

D+ = D[UD 3 UD 5 , D_=D 2 UD' 4 UD 6 . 

Since the fXj's and the $j's are well-behaved near A = and A = 00, respectively, we 
define M in terms of the ju/s in the regions D[, D2, -D3, and D' 4 , and in terms of $2 
and $3 in the regions D 5 and D 6 . The methodology of Fokas ( 19971 suggests making the 
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following ansatz for M: 



M(x,t,X) = < 



g[x,t) 
g[x,t) 
g(x,t) [ [j>3 



M2" ,,(12) 
a(A)'^3 

*4" „( 12 ) 

d(A) ' ^3 



(34) j£ 
d(\) 

a(X) 



- <I> 



(5) 



(6) ^ 



3(A) 



Ae As, 
Ae4 



(4.3) 



where 



d(A) = a(A)A(A) + b(X)B(X), 



(4.4) 



The definition of M in D[ Uf 2 U D3 U D' 4 , which involves the /i/s, includes the prefactor 
g(x,t). This prefactor is suggested by the relationship (2.8 1 between eigenfunctions of 



the Lax pairs (2.15) and (2.24), and its incl usio n implies that there exists a jump matrix 
J such that M_|_ and M_ are related as in (4.1). We introduce the following notation: 



J(x, t, A) = < 



r Ju 






A s 5j n L> 2 , 


j 2 = 


^3^4 




A e l» 2 n z> 3 , 








A G D 3 n L> 4 , 


J4, 






A G D' 4 nD[, 








A G L> 4 n L» 5 , 








A G L> 5 n D 6 , 




J4J5" 


■^6i 


\eD[n D 6 . 



(4.5) 



The jump matrices {J n }i can be determined from the various relations between the 
eigenfunctions. Indeed, algebraic manipulation of the equations (3.1 ) leads to expressions 
for the jump matrices {J n }t m term s of the spectral functions s(A) and S'(A). Similarly, 
algebraic manipulation of equation (3.6) leads to an expression for the jump matrix J$ 
in terms of the spectral function Soo(A). To find an expression for the jump matrix J5, 
we note that the relations (2.8), ( 2.14[ ), and (2.23) imply that two solutions fj, and <3? of 
(2.15) and (2.24), respectively, satisfy a relation of the form 



g(x,t)iJ,(x,t,X) = ^>{x 1 t,X)e' t6 ^ x ' t ^C{X)e l9{x ^ x) ' J \ 



(4.6) 



where C(A) is a 2 x 2-matrix independent of (x,t) and the functions 9(x,t, A) and 
oo (x, t, X) are defined by 



(x,t,X) 



x 

8A 



(A 



1 



>(x, t, A) = Xp(x, t) + 



t 



(4.7) 



In the particular case of fj, = [13 and $ = $2, equation ( |4.6| ) holds with C(X) — g(0, 0)s(A), 
g(x,t)fi 3 (x,t,X) = $ 2 (x,t,A) e -^( 2: ^ A )^ 3 (0,0) S (A)e l9 ( a; ^ A )^. (4.8) 



i.e 



Equation (4.8) together with (3.1a) and (3.6) provide the required relations between the 



/Ltj's and the $j's needed for determining J 5 . In summary, we arrive at the following 



11 




Figure 4 The image in the (y ,77) -plane of the half-line domain Q — {0 < x < 00, < 
t < T} under the map (4-- 11). 



expressions for the J„'s: 



where 



7° - 



Ji =e- ie&s J?, J 3 = e- t9a3 J? 



3 > 



7 5 = e - i9 ~ CT3 J 5 °e ie<T3 , 



1 

r(A) 1 



7° - 

J A — 




.6(A) 

a (A) a(A)a(A) 



W(AM0,0).(A)) ia fel(A ff 0))l ^ 

" ~ 1 (g(0,0) S (A)) 21 9(0,0)2 

aoo(A) 



Qoo(A)a(A) 



and r(A) is defined by 



r(A) 



a(A)d(A) : 



J 4 = e-^J 4 , 



7 U 



7° - 
j 6 — 



A G D 2 



_ (1 r(A) 

1 

1 



fcoc(A) 

a 00(A) 

1 



(4.9) 



boo (A) 

100(A) a = o(A)a 00 (A) / 



(4.10) 



4.2 RH problem for M(y, r], X) 

In the previous subsection we formulated a RH problem for M(x, t, A) in the Riemann 
sphere of the spectral parameter A. However, as noted above, this RH problem does not 
provide the solution of our initial-boundary value problem. Indeed, the jump matrices 
{Jn}\ involve 9^. The occurence of the function p(x,t) in 9^ implies that the RH 
problem cannot be formulated in terms of the initial and boundary data alone. To 
overcome this problem we make two important changes in the formulation of the RH 
problem: (a) We modify the jump matrix by adding appropriate exponential factors to 
the definition (4.3); (b) We introduce new variables (y, rj) by 



[x,t) (y,r)), y = p(x,t), V=n(x-t), 



(4.11) 



where p(x,t) was defined in (2.6 1. The jump matrix of the modified RH problem is 
explicitly given in terms of (y, 77, A) and can thus be formulated in terms of the initial 
and boundary data alone. 
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Define the 2 x 2-matrix valued function M(y, 77, A) by 



,/ I /, (12) 1 P iHv+t)c 

a I ^ // (12) ^ P* A (f+*) CT 3 



M(y,r),X) = < 



5 M 3 



(34) ^ \ iA(y+t)<7 a 



d(A) 



5 



(34) a4_ \ „iA(y+t) CT3 



_2 (f, 

a 00 (A)' V 3 



a(A) / 
( 5 M 



(6) *!> 



3 ) 



Qoo(A) 



XeD 2 , 

XtD 3 , 

XeD' 4 , 
Ae4 



(4.12) 



The map (4.11) is a bijection from Q — {0 < x < 00, < t < T} to a subset of the 
(y, ?7)-plane delimited by the image of dil under (4.11), see Figure |4j The image of dfl 
consists of three pieces given by 



(y,v)=2\l (1 - cos(u(x', 0))) V m(x', 0)dx', x 



< t < T 



(y, n) = \ (J (i + cos( u (o, t')))y/m(p,t>)dt', -t 

(V,V) = \ (p(0,T) + J (1 - cos(u(x',T))Wm(x',T)dx',x-T\ 
In particular, the map (4.11) is everywhere nonsingular: 





dy 




dx 

<1ll 






dx 


di 





|(1 — cosu)y/m ~(1 + cost^-^/m 



x>0 



= — -yrn < 0. 



(4.13) 
(4.14) 
(4.15) 



The expressions on the right-hand side of (4.12 1 should be understood as being evaluated 
at the point (x,t) corresponding to (y, 77) under (4.11). 

The form of the exponential factors e lX ^ v+t ^' y3 and e l ^ CT3 on the right-hand side of 
(4.12) is motivated by the fact that these exponential factors are analytic near A = and 

(4.16) 



A = 00, respectively, and by the relations 

e i\(y+t)<r 3 e i6o 3 _ e i(Aj/+ 3^)0-3 



The relations (4.16) imply that M satisfies the jump conditon 



M-(y,r],X) = M + (x,t,\)J{y,ri,\), A e D + n D 

M= {M + , XeD + , 
[M_, AeL>_, 

where the jump matrix J(y, i], A) is given by 

J„(y,77,A) = e-^ + ^J°(A), n = l 

( Ji Ae£>inz)2, 



.7: 



(4.17) 
(4.18) 

(4.19) 



J= < 



J 2 — J3J4 J\ 
J7 = J4J5 J& 



XeD 2 nD 3 , 
A e D 3 n D' 4 , 
A e D'^nD'^ 
A g D' 4 nD 5 , 

XeD 5 DD 6 , 

X e D[nD 6 . 
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ReX 



Figure 5 The contour for the Riemann-Hilbert problem in the complex X-plane. 



We can now prove the following theorem. 

Theorem 4.1 Let T < oo and let u (x), u\{x), < x < oo, and go(t), gi(t), < t < T, 



be given functions. Suppose that there exists a solution u{x, t) of equation (1.1) in the 
domain {0 < x < oo, < t < T}, which satisfies the initial conditions 

u(x, 0) = uq(x), Ut(x, 0) = u\(x), < x < oo, 

and the boundary conditions 

u(0,t) =g (t), u x (0,t) =gi(t), < t < T. 

Moreover, suppose that cos(u(x,t)) — 1 has sufficient smoothness and decay as x — > oo 
and that cos(go(t)) ^ —1 for < t < Then u{x 7 t) can be reconstructed from the 
initial and boundary values {uQ{x),u\{x) 1 gi){t) 1 gi{t)} as follows. 
Define the functions m(x, 0) and m(0, t) by 



m(x,0) = 1 + (u ax (x) + ui(x)) 2 , 
Define the functions p(x, 0) and p(0, t) by 
1 



m(0,t) = I + (gi{t) + g ot {t)f 



p(x,0) 



(1 — cos(uo(x')))\/m(x' , 0)dx' 1 



1 



^(0,*) = ^ / (l + cos( 5o (i')))vMO,i'X 



< x < oo, 
< t < T. 



Define /^(x, 0, A) and $3(1, 0, A) m terms ofuo(x), ui(x), and p(x,0) via the Volterra 
linear integral equations (2. Iffy and \2.21 ) evaluated att = 0. Define the spectral functions 
a(\), 6(A), &oo(A) 7 b co (X) by equations (3.2) and (3.1). Similarly, define ^1(2;, 0, A) in 
terms ofgo(t) and g\{t) via the Volterra linear integral equation (2.18) evaluated atx = 0. 



Define the spectral functions A(X), B(X) by equation (3.2) 
Let R > be such that 



R 2 > sup 



t-f 



o<t',t<T 4 J t ,(l + cos(5q(t)) V«(0, r)dr 



and define the sets {D n }\ by (2.20) and (2.32). Define D[ and D' 4 by (34). Let 
D+ = D[UD 3 UD 5 , D_ = D 2 U D' 4 U D 6 . 



1 If T = 00, we also assume that cos(u(0,t)) — 1 has decay as t — > 00. 
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Define d(X) and T(A) by (4-4) and (4-10). Assume that the possible zeros {kj}i of a(X) 
and {Xj}i of d(X) satisfy 



• a(k) has N simple zeros {kj}j =1 such that kj <E D[, j — l,...,ni, and kj 6 D-i, 
j = ni + l,...,N. 

• d(k) has A simple zeros {Xj}^ , such that Xj G Di, j = 1, . . . , A. 



None of the zeros of a(k) coincides with a zero of d{k). 



Then the solution u(x,t) of equation (1.1) is given in parametric form in terms of 
two real parameters y, rj such that 

T}>--; y>p{2r l7 0) for 77 > 0; y > p(0, 2|r?|) for - - < 77 < 0, (4.20) 



by 

where £ (y, 77) is defined by 



' 2lm (lp(2 V ,o) a (y^ r lW) +«o(2t?), 77 >0, 
2lm {lp(o,2\ n \) a (y'^)dy') + 5o(2|r?|), -f <r?<0, 



(4.21a) 



\f*o,m) V / l-4Im(a(i/',77)) 2 rf y ' + -f < 77 < 0, 

the function a(y, 77) is the unique solution of the Ricatti equation 

a y = a 2 -Ai\ lim XM 12 (y, 77, A) ) a - -, (4.21c) 

\ A— >oo / 4 

\ImA>0 / 

together with the initial conditions 

1 -iarcsinf ..o,(2^+. 1 (2^) 

a 0(277,0), 77) = e V ; rj>0, (4.21d) 

1 -iarcsinf , 91 ( 2 I 'il )+gQt ( 2 H) ^ T 

a(p(0,2|»7|),J7) = --e V V^TuTIiHR^hmV , -- < r/ < 0, 

and M(y, 77, A) is the unique solution of the following 2 x 2-matrix RH problem: 

• M is meromorphic away from the contour D + n D _ . 

• XTie possible poles of the first column of M occur at X = kj, j = 1, . .. , n\, and 
X = Xj, j = 1,...,A. The possible poles of the second column of M occur at 
X = kj, j = 1, . . . , n\, and A = Xj, j = 1, . . . , A. The associated residues satisfy the 
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following relations: 
Res[M(t/,7j,A)]i = 



a(kj)b(kj) 



(4.22a) 



Res[M(y,r),X)} 2 



-1 
a(kj)b(kj 



e-^+^^Miy^h, i = l,...,nx, 



(4.22b) 

Res[M(y,7?, A)]i =Resr(fc) e 2l(A ^ + ^ )<T3 [M(y, 77, A,)] 2 , j = l,...,A, (4.22c) 



Res[M(y,7 ? ,A)] 2 

A, 



Reside 2i{hv+ ^ )a3 [M{y, V ,\i)]u 

A A 



j = l,..., A. 

(4.22d) 



• M satisfies the jump condition 

M_ (y, r), A) = M + (y, 77, X)J(y, 77, A), A e Z?+ n ZA_, 
where M is M_ /or A G D_, M is A/ + /or A G 13+ , and J is defined by equations 



U-yj and (4.1ty . 
• M(y,r),X) =1 + 0(A), A 0. 
. M(y,r ? ,A) = / + 0(i), A -> 00. 

Proof. In the case when a(A) and d(A) have no zeros, the unique solvability is a conse- 
quence of the existence of a vanishing lemma. If a(A) and e?(A) have zeros, the singular 
RH problem can be mapped to a regular one coupled with a system of algebraic equations, 
see 



Fokas and Its (1996). 



(2008) implies that M as defined in (4.3) satisfies the residue condition 

1 



The residue conditions (4.22 ) can be p roved as follows. The general approach of Fokas 

3) satisfi 

e 2 ^ x ^[M(x,t,k 3 )] 2 , 3 = l 



R^[M{x,t,X)] 1 = w —- 



,n 1: 



where 9(x,t,X) is given by (4.7). The definition (4.12) of M and the relations (4.16) 
imply that M satisfies (4.22a). The other residue conditions follow similarly from the 
corresponding residue conditions for M. 



In order to prove (4.21), we note that the change of variables (4.11) implies that 

1 



d x = -(1 — cosu)^frndy + ^d v , d t = ^(1 + cosu)^/md. 



2^V 



In particular, 
Moreover, because 
we find that 



d x + d t = \pmdy. 

m = 1 + (u x + u t ) 2 = 1 + run 2 , 
1 



(4.23) 



l-ul 
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Thus, addition of the two equations in (2.241 together with (4.11) yields 



S,, + *A[0- 3| *] = ^ ( Uy + 



(4.24) 



Consider the particular solution of this equation given by 



$ = 



(5) 



MA)'* 3 



(5) 



Equation (|4.24|) implies that $ admits an expansion of the form 

1 



*(x,t,\)=I+ ^ + ^r^+°\ ^ 



A^oo, A e D 5 , (4.25) 



where < &^ 1 -' (x, t) and ^ 2 ^(x,t) are independent of A. Substituting this expansion into 
(4.24) we find by considering the terms of 0(1) that 



Thus, by construction of the RH problem, 



4 lim XMi 2 (y, rj, A) = u y + 

A— >oo 
Im A>0 



l-ul 



In view of the inquality 



T < Uy 



u x + u t 



< 1, 



y/l + (U x +U t )' 

we may define two functions Q(y,rj) and a(y,rj) by 

Q = -u- arcsin(u„), a = _^ e -"rc Si „(« H )_ 
These definitions imply that 



(4.26) 



(4.27) 



2Im(a) 



and that a satisfies the Ricatti equation 



1 - ul 



iQyOi 



(4.28) 



(4.29) 



Equation ( |4.21c| fo llows from ( |4.26[ ), ( |4.28| and ( |4.29[ ). 

Equation (4.25) implies that M = I + 0(1/ X) as A — > oo in D$. A similar argument 
shows that M = I + 0(1/ A) as A — > oo also in D e . 

In order to prove that a satisfies the initial conditions ( 4.21d| ) we note that the initial 
half-line {(x,t) = (2s,0)|s > 0} is mapped by (|4.11[) to the set {(y,rj) = (p(2s, 0), s)\s > 



0}. Together with the definition (4.27) of a this leads to 



a(p(2r),Q),r)) 



re <■ v '\x=2r,, 

'- t=0 



r?>0, 
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which yields (4.21d| for r\ > 0. The proof when —T/2 < r\ < is similar. 



Finally, equations (4.21a) and (4.21b) can be derived by considering the map (4.11) 



as the composition of the change of variables (1.2) with the map 



I'd 



me2£' +v,\v\ ~v),V 



(4.30) 



Thus, the variables {y,rf) £ Im(f2) satisfy (4.20). Using that the map (4.30) admits the 
inverse 



(y,ii) ^ i€,v) 



m l l 2 dy' + \r)\,r) 



together with the expression (4.28) for u y in terms of Ima, we find (4.21a) and (4.21b). 



□ 
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